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Abstract 



> 

■^j- ■ Recent calculations indicate that in the outer parts of neutron stars nuclei are rod- 

■ like or slab-like, rather than roughly spherical. We consider the elastic properties of 

these phases, and argue that they behave as liquid crystals, rather than rigid solids. 

We estimate elastic constants and discuss implications of our results for neutron 
00 . star behavior. 
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In the standard picture [1], the outer part of a neutron star is a solid from 
the bottom of the ocean of melted iron at a density p ~ 10 6 — 10 8 g cm" 3 
(depending on temperature) down to the boundary with the core of the star 
at a density of order the saturation density of nuclear matter, p s ~ 3 x 10 14 
g cm -3 , and matter consists of roughly spherical atomic nuclei arranged on a 
regular bcc lattice. Charge neutrality is provided by a background of electrons, 
and, at densities in excess of about 4 x 10 11 g cm -3 , there is a neutron liquid 
between the nuclei. However, recent work indicates that, in roughly half of the 
matter by mass in what is traditionally referred to as the "crust" of a neutron 
star, the nuclei are very aspherical, and have the form of essentially infinitely 
long rods or infinitely extended slabs (for a review, see Ref. [2]). Following 
Ref. [3], we shall refer to these phases as mesomorphic phases, since they have 
properties intermediate between liquids and solids. 

The possibility of nuclei in matter at sub-nuclear densities adopting very non- 
spherical shapes was first pointed out in the context of stellar collapse by 
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Ravenhall et al. [4], and confirmed by other authors using other models [5-8]. 
With increasing density, nuclei become first rods, then slabs. This is followed 
by two "inside-out" phases in which there are cylindrical or spherical "bubbles" 
in the nuclear matter. Subsequently the transition to the uniform liquid phase 
occurs. For neutron star matter with dripped neutrons, Lorenz et al. [9] found 
a similar sequence of phase transitions, and these conclusions were confirmed 
by other investigators [10,11]. 

Mechanical properties of the outer parts of neutron stars are important in 
attempting to understand many aspects of neutron star behavior, including 
glitches in the pulse repetition rate, magnetic field evolution, and neutron star 
models for 7-ray bursts [12-14]. Our purpose in this Letter is to discuss the 
elastic properties of the mesomorphic phases and estimate elastic constants. 

We shall concentrate on the phases with rod-like and slab- like nuclei, since they 
are expected to be the dominant ones. According to Ref . [9] , these "spaghetti" 
and "lasagna" phases constitute, respectively, roughly | and \ of the mass 
lying between the core and the ordinary solid crust, where nuclei are roughly 
spherical. In Ref. [11] these proportions are found to be about \ and |. Thus 
the fraction of mass in the bubble phases is relatively small. 

As a model, we shall assume uniform rods and slabs. For symmetric nuclear 
matter, Thomas-Fermi calculations [7] have confirmed that these simple con- 
figurations are the thermodynamically favorable ones. Clearly there is no in- 
crease in energy if rods or slabs are displaced in directions that lie in the plane 
of the slabs, or along the rods. Consequently there is no restoring force for 
certain sorts of distortion, and they thus have the elastic properties of liquid 
crystals. In the "lasagna" and "spaghetti" phases, microscopic calculations 
indicate that there is positional order in one and two directions, respectively, 
maintained by the Coulomb repulsion of rods and slabs [4-9] . Accordingly, they 
conform to the definitions of columnar phases and smectics A [3]. More com- 
plex positionally ordered phases (e.g., smectics C, cholesterics) are precluded 
by the symmetry of the equilibrium shapes of the nuclei. At the temperatures 
of neutron star interiors (~ 10 8 K) positionally disordered (nematic) phases are 
unlikely, since one would expect the ordering temperature of rods and plates to 
be comparable with the melting temperature for matter with spherical nuclei, 
~ 10 10 K. We emphasize, however, that the physical reasons for the spatial 
structure of the mesomorphic phases in the laboratory and in neutron stars are 
very different. For laboratory liquid crystals, the non-spherical shape of the 
molecules drives the tendency to form rod-like and slab-like structures, while 
in neutron stars, it is a spontaneous symmetry breaking brought about by 
the competition between the nuclear surface energy and the Coulomb energy. 
In the neutron star case, the basic objects, nucleons, from which structures 
are constructed are spherical, while in laboratory liquid crystals, the basic 
ingredients are non-spherical molecules. 
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To calculate energies of these phases we adopt a generalized liquid drop model, 
with bulk, surface and Coulomb energies. In the deformations that we study 
in this paper we assume that the total density remains constant. Distortion 
can lead to a redistribution of neutrons between nuclear matter and neutron 
matter, but this is small because bulk energy densities are large compared 
with those of surface and Coulomb energies. Therefore we may assume that 
the fractions of the total volume occupied by nuclear matter and neutron 
matter, and their local densities remain constant, and consequently only the 
Coulomb and surface energies are altered. 

We begin by estimating the elastic constants for the layered phase. Since there 
is complete rotational symmetry about the axis perpendicular to the layers, 
which we denote by Oz, this phase is similar to a smectic A liquid crystal, 
and the energy density due to deformation may be written in the form [3] 



dU 1 ,n 



+ : y(Vi«) 2 , (1) 



where u is the displacement of layers in the z direction. The first term is 
associated with a change of interlayer distance (Fig. la). The shear shown in 
Fig. lb is equivalent, due to rotational invariance, to reducing the interlayer 
spacing, an effect taken into account by the second term in the square brackets. 
The second-order elastic constant K\ is associated with splay deformations. 
As a first case, we consider a distortion that maintains the planar character of 
the slabs but changes the layer spacing, which we denote by 2r c . The surface 
energy per unit volume, E a scales as r^ 1 , while the Coulomb energy per unit 
volume, Eq scales as r 2 , and therefore in equilibrium these energies satisfy 
the condition E a Q = 2Eco [2,4], which determines r c . Here the subscript "0" 
denotes equilibrium values. From the same scaling law we conclude that a 
small departure from equilibrium associated with a change 5r c of r c results in 
an increase of energy 

5E a+c = 3E C0 (5r c /r c ) 2 . (2) 



Applying Eqs. (1) to the same deformation, we obtain Ed = \B(5r c /r c ) 2 , 
whence it follows that 

B = 6E C0 - (3) 



In order to estimate K\, let us consider a harmonic perturbation u(x) = 
Uocos(kx), assuming that its wavelength is large: (kr c ) 2 <C 1 and (kuo) 2 <C 1. 
Then the deformation energy averaged over x is 

(E d ) = (ku o r + \K x k'ul (4) 
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Fig. 1. Linear deformations. Cross sections of the slabs and rods in equilibrium are 
indicated by full lines, and, after deformation, by the hatched areas, a: compression 
perpendicular to the slabs (rods); b: shear; c: transverse shear in the columnar phase; 
the primitive translation vectors (ai, a2) and their counterparts in the reciprocal 
lattice (ki, k2) are also shown. 

The first term in this expression dominates when uq ^> A, where A = {K\j B) 1 / 2 . 
To order k A , the relative change in the surface energy averaged over x is 



(AK> 



E, 



(ku f 



0-0 



64 



(5) 



If one includes curvature contributions to the surface energy one finds that the 
first order terms average to zero, while the second order ones give an additional 
contribution bk 4 u 2 ] /2 to this equation, where b is the ratio of the second-order 
curvature coefficient to the surface tension. For symmetric nuclear matter, 
according to the estimate of Bennett and Ravenhall [15], b —(0.4 fm) 2 . 

It is convenient to express the charge density of nuclei as a Fourier series, and 
for the distortion given above one finds 

p(x,z) =J2p mn exp( y im7i-z/r c + inkx). (6) 

mn 



Here p mn = p p (mn)' 1 sin(m7rrjv/r c )(— i) n J n (miruo/r c ), p p being the charge 
density of protons inside nuclei, r N the halfwidth of the slab, and J n the 
Bessel function of order n, and the summation is over all integers. The spatial 
average of the Coulomb energy is 

{EC) = £ Wf + {rnn/Z? + (nkr (?) 



where kxF is the Thomas-Fermi screening wavenumber, and the prime on the 
summation indicates that the term with m = n = is excluded. In the absence 
of distortion (-u — 0), the screening correction can be neglected in evaluating 
the Coulomb energy, since k\ F r^ 2 , as argued in Refs. [7,16], and the 
summation over m in Eq. (7) yields Eco = 27r (p p rN) 2 {l — w) 2 /3, in agreement 
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with Ref. [4]. Here w is the fraction of space occupied by nuclear matter, and 
w = rjv/r c for the smectic case. From the equilibrium condition E a0 = 2Eqo 
and the fact that the surface energy is given by E a0 = wcr/r N , where a is the 
surface tension of the interface between nuclear matter and neutron matter, 
it follows that in equilibrium r N = (3wo-) 1/3 (4irp 2 p (l - w) 2 )~ 1/3 . Thus the 
Coulomb energy in equilibrium varies as the two thirds power of the surface 
tension. 

Although k,TF may be large compared with k, the effects of screening may 
be neglected also in calculating the change in the Coulomb energy due to 
distortion, (AE C ) = (E c ) — E C q, because p Qn = for any nonzero value of n. 
Thus we find 

(AE C ) = { Pp r N f £(-1) (kr c /n)V, (8) 

where 

%(0 = ^E m 4 +2j E ^ Jn(mO, (9) 

m=l n=l 

with ^ = nuo/r c . In particular, we obtain ai(£) = 7r^ 2 (1 — w) 2 /3 and a2^) = 
it£ 2 (1 - w) 2 [(1 + 2w - 2w 2 ) tt 2 /45 + i 2 /A}. The leading Coulomb contribu- 
tion to the energy due to distortion is — E C o(ku ) 2 /2. This exactly cancels the 
corresponding surface term in Eq. (5), as it must on general grounds. By com- 
paring the sum of the Coulomb (Eq. (8)) and surface (Eq. (5)) contributions 
to the energy due to distortion with the general expression (4) we find 

K 1 = 2E C0 (l + 2w -2w 2 )r 2 J15. (10) 

The curvature correction to this expression is 4bEco, which may be neglected 
since r c is typically ~ 10 fm. 

In the calculations above we assumed that the thickness of a slab measured 
along Oz remains unperturbed. In reality, the slab will be thinner near the 
extrema of u(x) and thicker near the extrema of du/dx. Allowance for this 
effect, however, yields a contribution to (E d ) containing an additional small 
factor ~ (k 2 + k^ F ) r 2 relative to the terms already considered. 

Finally, we remark that, even though at finite temperatures strict long range 
positional order of layered phases will be destroyed according to the standard 
Landau- Peierls argument, our theoretical estimates of elastic constants should 
be a good approximation at low temperatures. 

Now we turn to the columnar phase. We assume the rods to have a circular 
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cross section of radius and define r c by the condition that the density of 
rods per unit area is 1/irr 2 . Then the scaling laws that lead to Eq. (2) hold 
also for this case. 

The lowest energy configuration is one with rods on a two-dimensional trian- 
gular lattice. We choose Oz to be the axis of the symmetry. Displacements 
are then described by a two-dimensional vector u = (u x ,u y ), and the energy 
of deformation is 




The first two lines of this formula reproduce Eq. (7.28) of Ref. [3]. The con- 
stant B is associated with uniform transverse compression or dilation, C with 
transverse shear (Fig. lc), and K 3 with bending. 

The third line contains higher-order terms. The last one may be important 
when amplitudes of longitudinal deformations exceed (K 3 / B") 1 / 2 . The term 
with B' provides a non-linear coupling of transverse and longitudinal defor- 
mations, which causes, among other effects, a breakdown of linear elasticity 
and hydrodynamics at large scales (see §§8.3-8.4 of Ref. [3]). In the case of 
smectics, an analogous coupling is provided by the cross term coming from 
the expression in square brackets in Eq. (1). 

The terms B, B', and B" stem from the energy increase due to the change of 
the cross section of the unit cell. A longitudinal shear du x /dz = A accompa- 
nied by uniform transverse compression du x /dx = du y /dy = A'/2 maintains 
the triangular lattice but leads to a change of the lattice spacing propor- 
tional to (1 + A') 1/2 /(l + A 2 ) 1 / 4 . Then Eqs. (2) and (11) yield B = \E C o, 
B' = — f-Eco, and B" = |i?co- 

The two other constants are associated with deformations that are not a simple 
scaling of the triangular lattice. We estimate them in a way analogous to that 
for smectics. Imposing a perturbation of the form u(z) = u cos kz and Fourier 
transforming, we arrive at the expression 



(F _ ^/ 2%wp 2 p J 2 (k lm ■ u ) / 2J 1 (k lm r N ) \ 
L k 2 TF + k 2 m + (nk) 2 { k lm r N ) 



(12) 



6 



where k; m = /ki + mk 2 is a reciprocal lattice vector. In the absence of trans- 
verse deformation, kx and k 2 are vectors of length (87r/v / 3) 1/ ' 2 ^ ( r 1 , with an 
angle 7r/3 between them (Fig. lc; we shall assume that k 2 is directed along 
Oy). With uq = 0, only the term n = survives, and neglecting k TF , we re- 
cover Eq. (10) of Ref. [6] for Eco- The analytic formula obtained by replacing 
the hexagonal unit cell by an equivalent cylinder [4], 

E co = (tt/2) ( Pp r N ) 2 w (Hl/w) -l + w), (13) 

turns out to be very accurate at the values of w at which the "spaghetti" phase 
is expected^]: it underestimates E C o for the hexagonal cell by less than 0.6% 
at w < 0.35, which is several times smaller than the difference due to inclusion 
of a realistic screening wavenumber kTF ~ 0.4 r~ l . Since for the columnar case 
the surface energy is E a = 2wa/riy, one finds r^v = (2a) 1 / 3 (7rpp[ln(l/ 'w) — 1 + 
w])~V 3 m equilibrium. 

A transverse shear du x /dy = A changes k± y by (27TA/3) 1 / 2 A/r c (Fig. lc). Nu- 
merical summation of the series (12) and identification of the term quadratic 
in A yields the elastic constant C, shown in Fig. 2a. In the relevant range 
of w, it can be approximated as log 10 (C / Eco) ~ 2.1 (u> — 0.3) (dashed line). 
(This estimate of C is only a first approximation, since we have kept the cross 
sections of rods circular, whereas in reality they can adjust their form to the 
shear) . 




h Ac) 2 (*■„ Ac) 2 

Fig. 2. Shear constant C and bend constant K3 of the columnar phase. 

Finally, using a finite uo in Eq. (12), performing the summation and identifying 
the term proportional to k 4 u^, we determine the bend constant K 3 , which is 

1 The columnar phase is favored for w between a lower limit of 0.15-0.20 and an 
upper limit of 0.30-0.35 (and an upper limit of w for the layered phase is about 
0.60-0.65) [4-7]. 
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plotted in Fig. 2b. In the physically relevant range of w, K% = (0.064 — 
0.067) Eco^c- The term proportional to /c 4 Mq yields for B" the value obtained 
earlier. 

The fact that matter in neutron stars near the boundary of the core has the 
elastic properties of liquid crystals rather than a crystalline solid will have 
important consequences for a number of aspects of neutron star behavior. 
First, the maximum elastic energy that can be stored in the crust will be 
reduced, and this should be taken into account in models of glitch phenomena 
and starquakes that depend on the deviations of the figure of the star from that 
for a fluid (see, for example, Ref. [12]). It is not possible to make a quantitative 
estimate of the reduction of the elastic rigidity of the outer parts of the star 
without a detailed model for the orientation of the liquid-crystal phases, but 
since these phases are estimated to make up roughly half of the matter by mass 
outside the core, one might expect the effective elastic rigidity to be reduced 
by roughly a factor of 2. Since the "plate tectonics" of liquid crystals is likely 
to be very different from that of crystalline solids, a second problem where 
the elastic properties will be crucial is in models of the evolution of neutron 
star magnetic fields that invoke such processes [13]. Other applications are 
to the propagation of elastic distortions, and to the energy of defects in the 
liquid-crystal structure. 

In our discussion above we assumed that in the absence of deformation the 
spacings of rods and slabs had their equilibrium values. However, in reality, this 
may not be the case, since to alter the spacing requires a major rearrangement 
of the proton distribution. Just how easy it is to do this depends on how 
slabs and rods are connected to each other, and the related question of what 
defects are present. Should the spacing be smaller than its equilibrium value, 
there will be a positive contribution to the elastic energy of the form {E a /2 — 
-E'c)(V ±u) 2 , which gives a restoring force linear in the deformation. Should the 
spacing be larger, this contribution will be negative, which gives a tendency 
to rotation and leads to spontaneous deformation analogous to the Helfrich 
effect in conventional liquid crystals [3]. Eventually the spacing will approach 
equilibrium, either by spontaneous reconnection of rods and slabs or by motion 
of defects. 

To summarize, the above considerations indicate that matter near the bound- 
ary with the core of a neutron star has the elastic properties of a liquid crystal, 
rather than a conventional solid. To distinguish this region from that further 
out in the star, we consider it appropriate to refer to it as the "mantle" rather 
than as part of the "crust". It is clear that the properties of neutron stars 
need to be reconsidered in light of this new understanding. 
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